Depletion theory and the precipitation of protein by polymer 
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The depletion theory of nanoparticles immersed in a semidilute polymer solution is reinterpreted 
in terms of depleted chains of polymer segments. Limitations and extensions of mean-field theory 
are discussed. An explicit expression for the interaction between two small spheres is derived. The 
depletion free energy for a particle of general shape is given in terms of the capacitance or effective 
Stokes radius. This affords a close to quantitative explanation for the effect of polymer on protein 
precipitation. 
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It is difficult to set up theories of the polymer distri- 
bution near surfaces because the computed polymer in- 
homogeneity is sensitive to the nature of the approxima- 
tions introduced. Pierre-Gilles de Gennes devoted con- 
siderable time and effort to trying to understand these 
problems starting with his survey [l| of forty years ago, 
which is still illuminating to read, and culminating in 
papers containing now classic ideas like self-similarity [2] 
and the proximal exponent [3|. His concise, powerful note 
[^ on polymer depletion by a small sphere was strangely 
neglected by the colloid community for a long time until 
the peculiarity of depletion on nanoscales was reassessed 
merely a decade ago [a, |g] . There has been a flurry of ac- 
tivity in the statistical physics of nanocoUoids immersed 
in polymer solutions in what has been termed the pro- 
tein limit (see e.g. 7, 8, 9, 10, 11] and references therein). 
Here, I would like to emphasize simple aspects of polymer 
depletion by nanoparticles in the spirit of ref . jl2| . 

Let us recall the argumentation introduced by de 
Gennes YS\ to compute the free energy of depletion in- 
volved in immersing a nanosphere into a semidilute poly- 
mer solution. The solvent is not just "good" but needs 
to be really "excellent" (see below), i.e. the excluded vol- 
ume P between the Kuhn segments equals A^ where A is 
the segment length. If the radius a of the sphere is larger 
than A, it is plausible to assume that a and the polymer 
correlation length ^ are the only relevant length scales in 
the problem. The latter is given by [l3| 
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where the concentration cq is the number of polymer seg- 
ments per unit volume. For a nanosphere dissolved in a 
semidilute solution of low concentration, one readily has 
a ^ ^. De Gennes then argues that there is a volume 
of order a^ ~ independent of ^ - surrounding the sphere 
from which polymer is depleted \^. Hence, the number 
of segments depleted is of the order of a^ cq- Since the 
free energy of depletion Fi must be proportional to this 
number, i.e. must be proportional to the concentration 
Co, one concludes that 3 



Fi 
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Here, fee is Boltzmann's constant and T is the tempera- 
ture. 

Here, I invoke a different type of scaling argument 
to derive Eq.(2) because this will allow a direct assess- 
ment of the other aspects of depletion in the analysis 
below. The number of polymer segments depleted from 
the vicinity of the sphere is o^ cq so one might be tempted 
to think naively that the free energy of depletion could 
be something like Fi ~ a^ cokBT. But this disregards 
entirely the fact that the segments are all connected 
(and are actually all on one single polymer chain because 
0^^). A depleted test segment is connected to h others 
where a~ h^'^ A, in view of the excluded-volume effect. 
Therefore, the number of degrees of freedom is reduced 
by a factor h 



Fi ~ a^ Co h-'^ kBT , 
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which agrees with Eq.(2). Note that this derivation is 
valid in the mean for the effective number of degrees of 
freedom is actually less than unity in Eq.(3). 

In the following I discuss several problems concerning 
the depiction interaction between nanospheres and poly- 
mer and its usefulness in explaining the precipitation of 
proteins by polymer. 

Water is often the solvent of choice in experiments 
concerning the thermodynamic properties of polymer- 
nanocoUoid or polymer-protein mixtures. Although 
water-soluble polymers dissolve readily in aqueous solu- 
tion, the solvent must often be regarded as "intermedi- 
ate" (/?<A3) rather than "excellent" [i|. The polymer 
chain then interacts with a nanoparticle in a quasi-ideal 
manner [12l |. 

In effect, the excluded-volume parameter z = h}/"^ P/A^ 
may remain smaller that unity if the nanoparticle is not 
too large. In that case, a string of depleted segments 
behaves like a Gaussian chain: h~a^/A'^ if the particle 
is a sphere. This would imply the condition 12:] a< A^//? 
which may be easily met in practice. Eq.(3) then leads 
to 111 
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1 
0.614 



1.1 
0.585 



1.2 
0.558 



1.3 
0.532 



1.5 
0.486 



2 
0.395 



4 
0.222 



6 
0.1538 



10 
0.0952 



TABLE I: Scaled depletion energy between two nanospheres as a function of their separation 



where fci is a numerical coefficient and ^id — 1/ {A^cq) is 
a quasi-ideal correlation length. The latter relates to a 
state in which the polymer chains are supposed hypothet- 
ically ideal. Within the same approximation, one may 
argue in favor of a self-consistent field picture at a ^ ^ 
which leads to a Laplace equation 12] for ^(r) where the 
inhomogeneous polymer segment density c(r) = co^(r)^ 



A*(r) =0. 



(5) 



The solution to Eq.(5) with ^ = at the boundary of the 
sphere is 



*(r) = 1 - 



(6) 



This leads to a value fci ~ 27r/3 for the coefficient in 
Eq.(4). 

Nevertheless, this SCF point of view cannot be entirely 
correct. At some distance r, from the sphere, excluded 
volume effects must come into play in Eq.(5) but then the 
argumentation for a mean-field approach also becomes 
weak. From renormalization theory [6| , we know that in 
the case 13 = A'^ we have asymptotically 



*M-i-(-) 



(7) 



with exponent x ~ 4/3. When the solvent is intermedi- 
ate {(i <^A^), we again adduce reasoning based on the 
parameter z above to show that Eq.(6) is only valid for 
r<r^=A'^ I (3. Beyond r», Eq.(6) must join smoothly to 



*(r) ~ 1 
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This reduces to Eq.(7) if a>r*. 

The SCF theory for the polymer distribution near a 
surface is incorrect for a number of reasons. First, the 
correlation length ^ is given by a wrong power law which 
de Gennes |3| proposed to amend by changing the expo- 
nent in the excluded volume term in the SCF equation. 
Nevertheless, the equation remains purely diffusive so it 
cannot mimic the segment distribution about a sphere 
expressed by Eq.(7). A further amendment could be to 
introduce a fractional SCF equation. In the vicinity of 
the sphere, this would reduce to 
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which is a stationary generalized diffusion equation of 
fractal order D and with a diffusion coefficient r~^. (For 
a discussion of fractional diffusion equations, see ref. 
[15|). Eq.(7) imposes the constraint D — Q = 10/3 on 



the exponents D and 0, which leaves one of them to be 
suitably chosen. However, a fractional SCF equation is 
still not entirely satisfactory. Density fluctuations are 
merely accounted for in a preaveraged sense. 

It is interesting to note that there is a relation between 
the depleted concentration about a sphere c{r) — co and 
the pair correlation function g{r) pertaining to a single 
chain in the bulk [I^ (see Eqs.(6) and (7)) 



c{r) - Co ~ -Co a^ h ^ g{r) . 



(10) 



This is independent of the strength of the excluded- 
volume effect. The structure of the depletion hole is 
analogous to the pair correlation structure of the polymer 
removed provided the latter is normalized by the effective 
number of degrees of freedom depleted by the sphere. 

It is of interest to consider the interaction between two 
nanospheres separated at distance R in the quasi-ideal 
limit. One needs to solve Eq.(5) with boundary condi- 
tions ^ = at their surfaces and \1/ = 1 at infinity. In 
the electrostatic analogy, the two spheres are grounded. 
This is not a trivial problem [10]. Using the method 
of images [l3|, one proceeds as follows. A positive test 
charge is placed at the center of the first sphere. The 
potential ^ on the second sphere is brought to zero by al- 
ternately adding appropriate image charges on the center 
line between the two spheres: these are negative within 
the second sphere but positive in the first. Next, a simi- 
lar test charge of positive sign is placed at the center of 
the second sphere. The potential on the first sphere is 
now rendered uniform by again adding image charges on 
the centerline but the signs are interchanged. The po- 
tential on both spheres is identical and uniform. Thus, 
we derive the capacitance of the two spheres as a series 
expansion which may be written in a less cumbersome 
manner via the method of difference equations. The de- 
pletion free energy turns out to be directly related to the 
capacitance of a particle [13] via Green's first identity 
which finally leads to 
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Eq.(ll) disagrees with an expression quoted by Hanke et 
al. without proof [l8|. The function 2 — C2 / {'in a) is very 
slowly varying (see Table 1) and has a maximum equal 
to 2 — 2 In 2 as the spheres touch. Here, the units for 
capacitance have been chosen in such a way that Ci = 
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FIG. 1: The slope of the hnear plot of ^"logS' versus w 
as a function of the effective Stokes radius of various pro- 
teins. The data are taken from ref. [26l |. The proteins are: 
lysozyme (Lys), a-lactalbumin (a-Lac), chymotrypsin (Chy), 
human serum albumin (HSA), human 7-globulin (IgG), al- 
dolase (Aid), thyroglobulin (Thy), human fibrinogen (Fib) 



47ra for a single sphere. Inevitably, Eq.(ll) breaks down 
at separations R beyond r* as the excluded- volume effect 
starts to play a role. Nevertheless, the decay of F2 with R 
is not fast enough for a finite computation of the second 
virial coefficient. A cut-off at i?~^ must be introduced 
as has been discussed by Eisenriegler [1^ in the excellent- 
solvent case (/3 = A^). The opposite limit of large spheres 
near ideal polymers has been addressed by Tuinier et al. 

How well does an expression like Eq.(4) work? An 
important phenomenon bioengineering is the precipita- 
tion of proteins by inert polymer. This has been studied 
for a long time [2l|, HI, ^ IH, Hi] but the most thor- 
ough quantitative study is that of Atha and Ingham [2^ . 
They determined the solubility S of a host of proteins 
as a function of the polyethylene glycol (PEG) added to 
the suspension. The concentrations of PEG were well 
into the semidilutc regime. The logarithm of the solubil- 
ity turns out to be a purely linear function of the PEG 
concentration and the resulting slopes are a monotone 
function of the protein radius, the latter being an equiv- 
alent quantity derived from the diffusion coefficient via 
the Stokes-Einstein relation (see Fig. 1). The proteins 
may deviate significantly from an ideally spherical shape 
so let us account for this fact. First, we know that for a 
compact particle of general shape, Eq.(4) may be gener- 
alized to 



in terms of the capacitance Ci [I4I. Next, Hubbard 
and Douglas have shown analytically and by simulation 
[27l . l28j that the Brownian friction coefficient of a par- 
ticle of general shape is directly proportional to its ca- 
pacitance to an excellent approximation. (It is also well 
to recall that the capacitance itself is essentially propor- 
tional the the particle's surface area [23,[30|.) Since the 
chemical potential of a protein given by the sum of Fi and 
k^T InS must be a constant, experiment should conform 
to the expression 



A 1" log 6* 
Aw 



= -0.051 as : 



(14) 



valid for PEG solutions where w is given in % weight 
per unit volume and as is the Stokes radius of the pro- 
tein in nm. (For data on PEG under theta conditions, 
see ref. [3l|.) In Fig. 1, we have plotted Eq.(14) but 
with an adjusted coefficient -0.043=- (0.23/0.27) x 0.051 
because PEG4000 isn't quite long enough to be char- 
acterized as infinitely long (see Fig. 3 and Table 2 of 
ref. [201). Although Eq.(14) overestimates the impact of 
polymer a bit, it is remarkably how well a simple SCF 
model works. At high protein radii, one expects signif- 
icant upward deviations from linearity owing to Eq.(2) 
but the trend is the reverse which is puzzling. Appre- 
ciable attractive forces between a protein and PEG may 
sometimes exist as has been suggested by Bloustine et al. 
[22] (see also ref. [331) ^'^^ could be the cause of these 
anomalies. 

In this note and in previous work [3J], we conclude 
that linear depletion laws for thermodynamic quantities 
are valid up to quite high protein concentrations. The 
linearity is in accord with scaling and SCF arguments 
for nanoparticles. As yet there is little evidence for 
interaction terms as given by Eq.(ll) for example. It 
is quite possible that repulsive forces between proteins 
(hard core, electrostatic) are compensated by attractive 
forces (adhesive, depletion) in such a way that quasi-ideal 
conditions apply. For a theory of this effect, see ref. [39 ]. 
The concentration of electrolyte added to the mixtures 
in precipitation experiments [26| would suggest that we 
are in such a regime. 
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